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In writing a covariant effective action for single field inflation, one is allowed to add a Gauss- 
Bonnet and axion-type curvature couplings. These couplings represent modifications of gravity, 
and are the unique higher-curvature terms that lead to second order equations of motion in four 
dimensions. In this paper we study the observational consequences of such couplings for models with 
large non-gaussianities. Our focus is on the Gauss-Bonnet term. In particular, we study an effective 
action where the scalar Lagrangian is a general function of the inflaton and its first derivative. We 
show that, for large non-gaussianities, one can write /jvl in terms of only three parameters. The 
shape of /ml is also studied, and we find that it is very similar to that of k-inflation. We show that 
the Gauss-Bonnet term enhances the production of gravitational waves, and allows a smaller speed 
of sound for scalar perturbations. This, in turn, can lead to larger non-gaussianities which can be 
constrained by observations. Using current WMAP limits on /ml and the tensor/scalar ratio, we 
put constraints on all parameters. As an example, we show that for DBI inflation, the Gauss-Bonnet 
coupling leads to an interesting observational window with both large /ml and a large amplitude 
of gravitational waves. Finally, we show that the Gauss-Bonnet coupling admits a de-Sitter phase 
with a relativistic dispersion relation for scalar perturbations. 

PACS numbers: 



I. INTRODUCTION 

Cosmology has entered an era of unprecedented progress. High precision measurements of the cosmological param- 
eters have led to a coherent picture of the history of our universe that seems to favor the inflationary paradigm 
Moreover, a future detection of large non-gaussianity in the cosmic microwave background (CMB) would falsify the 
simplest inflationary scenario, namely, single held slow-roll inflation [2|, y, |4[. 

On the theoretical side, there has been great activity in trying to produce large non-gaussianities in single and 
multiple-held inflationary models. For single field inflation, large non-gaussianities are easiest to produce in models 
with a small speed of sound (see e.g. [1, Hj]). On a parallel set of developments, there has been recent interest in 
developing a systematic effective field theory of single field inflation 0, H|. In ref. 0, such approach was applied 
directly to the Lagrangian describing the perturbations around the inflationary solution. The effective action can be 
viewed as an expansion in powers of (g m + 1) and the extrinsic curvature K a b of the constant time hypersurfaces. 
Such approach is quite general, and provides a straightforward way of calculating all CMB observables directly from 
the effective action for the fluctuations. 

On the other hand, one would like to understand how the various terms in the effective action for the fluctuations 
relate to the effective action of the inflaton itself. A method to build such an effective action was introduced by 
Weinberg in [8[. In this approach, one considers all marginal and irrelevant operators involving the inflaton and the 
metric. Among these terms, there are higher curvature invariants coupled to the inflaton. Generically, such terms will 
contain higher time derivatives on the helds which need to be eliminated using the hrst order equations of motion. 
Otherwise, one would be propagating more degrees of freedoms than intended. Weinberg showed that, after such 
eliminations and to leading order in the derivative expansion, the resulting action for the inflaton takes the familiar 
k-inflation type form plus two extra couplings between the inflaton and the Weyl tensor. Such extra couplings can 
be written instead in terms of the Gauss-Bonnet tensor and an axion-type coupling. In this way, one has an effective 
action that leads to second order equations of motion explicitly. 

One can then ask if it is possible to re-sum such an expansion. Moreover, one would like to write down a general 
local action for the inflaton coupled to gravity that leads to second order equations of motion for all fields. By abuse 
of notation we will call such re-summation a "UV completion" of the effective theory. An advantage of having such 
action is that one can then build directly the low energy effective action for the fluctuations as in [7j, have a clear 
physical interpretation of the various couplings, and a better assessment of their relative importance. If we also insist 
in preserving general covariance in the UV completion, it is not hard to see that an obvious candidate for such action 
is ^ 

^R + P(X,cf>) + V 1 (cb)E 4 + V 2 (<b)e abcd R ab ef R cdef . (1) 
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where E4 is the Gauss-Bonnet combination, 

and X = — iV a (/)V a is the kinetic term. In writing Q we assume that we work in the Einstein frame. The 
second term in the action is the familiar k-inflation type [l0( . The last two couplings in the effective action would be 
topological invariants in four dimensions if both potentials Vi, V2 were constant. This is the reason why they lead to 
second order equations of motion for general Vi. Note that if these potentials depended on X, one would end up with 
equations of motion depending on more than two time derivatives. 

The last two terms in the action |T]) represent modifications of Einstein's gravity, and as such, they have been 
studied numerous times (see e.g. [Hj]). Moreover, such couplings are known to arise in string theory The last 
term in (|TJ) is of the axion-type and it was studied long ago in [141 ] , where it was shown that such coupling does not 
affect the evolution of scalar fluctuations to quadratic order. We have verified that this is still true to cubic order. 
Therefore, we discard this coupling in what follows 25]. The Gauss-Bonnet term, on the other hand, has been studied 
many times in the context of Dark Energy (e.g. For other studies in the context of early cosmology see 

This term does contribute to scalar fluctuations and it is the main focus of this work. 

The purpose of this paper is to study the observational signatures of the Gauss-Bonnet coupling in the context of 
inflation. Moreover, our main interest will be in models will large non-gaussianitics. In the following, we compute 
the non-gaussianity parameter f^L using the action ([I]) in the limit of a small speed of sound. Moreover, we perform 
such calculation to leading order in the slow roll parameters, but to all orders in the "strength" of the Gauss-Bonnet 
coupling, defined as 

g = 8V{(<f>)H4, , (3) 

where H is the Hubble parameter. In making the calculations, we assume that the parameter g is slowly varying in 
time. We find that for large-nongaussianities, //vz, can be written in terms of only three parameters. We study the 
shape of of Jnl as a function of these parameters. We find that the shape of /jvl is always very close to that of 
k-inflation, even in the limit g — > 00. We discuss in which cases deformations from this shape might be observed. 

The spectrum of gravitational waves is also studied, and we find an enhancement due to the Gauss-Bonnet term. 
Using WMAP limits on the equilateral Jnl and the tensor/scalar ratio, we put constraints on the different parameters. 
For the particular case of DBI inflation, we show that one is left with only a two-parameter family of $nl- In this 
case, one can put a more precise constraint on the Gauss-Bonnet coupling: 

Sdbi < 3 . 

We also discuss implications of the Gauss-Bonnet coupling on the Lyth bound in the context of DBI inflation [l5| . 

An interesting aspect of the Gauss-Bonnet coupling is that scalar perturbations are non-trivial even in a de-Sitter 
background, just like the do in the Ghost Condensate (lr| . However, we find that in our case, quadratic scalar 
fluctuations have the familiar relativistic dispersion relation uj ~ k instead of the non-relativistic one u> ~ k 2 of the 
Ghost Condensate. This matches perfectly with the new de-Sitter limit found in [3] using the effective action of the 
fluctuations. Therefore, the Gauss-Bonnet coupling is precisely the modification of gravity that leads to such limit. 

The paper is organized as follows. In section II we study the equations of motion for the background solution and 
the quadratic scalar fluctuations. We also point out the different limits used in the calculations. In section III we 
calculate $nl and study its shape. In section IV we calculate the gravitational wave spectrum. We discuss the various 
constraints on the parameter space. In section V we study the case of DBI inflation. Finally, we close with some final 
comments and future directions in section VI. 



II. BACKGROUND SOLUTION AND QUADRATIC FLUCTUATIONS 

The equations of motion for the homogeneous background that follow from the action |T]) can be written as, 

E = 3H 2 (l+g) , (4) 
E = -3H(E + P) + 3H 3 (l-e)g , (5) 

where 

E = 2Xd x P-P, X = \i>\ e = ~Jp' 
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and g is defined in Eq. ([3]). 

Since we are interested in isolating the contribution to Jnl of the Gauss-Bonnet coupling, we will set e = 0. That 
is, we will work in a de-Sitter background [26j]. The existence of such limit will be established in the next section. 
Moreover, we will assume that the parameter g varies slowly with time. Therefore, to the first approximation we can 
consider it to be a constant. Note however, that g itself can be large. Therefore, we will do our calculations to all 
orders in g. 

Under these assumptions, one can easily show from (0| and © that, 

X8 X P = Ui 2 g . 

This will be a very useful relation in what follows. The speed of sound of this model is given by, 

Since we are interested in the limit of large non-gaussianities, we will be working with a small speed of sound. In this 
case, one can eliminate time derivatives of cj> by 

where we have defined 

In the limit of small speed of sound, we do not need to assume that 4>/H 3> 1, but only that <f>/H is much larger 
that any of the slow roll parameters. By slow roll parameters we mean any of the parameters that encode the time 
evolution of the solution such as e, r/ = e/ He, etc. 

In this case, one can show from the equations of motion ([4]) and (0 that derivatives with respect to the field <j> are 
suppressed, e.g. 

d^P (slow roll) 

If WT~ <<: ' 

Similar limits can be shown for other quantities involving derivatives of the scalar field like, e.g., d^dxP- 
Let us now define the parameters 

Note that £ is a slow roll parameter, and by assumption, it must be small. It is then easy to show that, 

Therefore, in the slow roll limit, rji is small compared to g. This means that we will be able to ignore higher derivatives 
of the potential V{(<j>) in the following. 

Note that the slow roll conditions are necessary to ensure an almost scale-invariant spectrum of scalar fluctuations. 
Nevertheless, we will see that the value of g does not affect the scalar tilt, and hence it can be much larger than any 
of the slow roll parameters. 



A. Quadratic Fluctuations 



To derive the action for the fluctuations we will use the ADM formalism, where we write the metric as 

ds 2 = -N 2 dt 2 + h ij (dx i + N l dt){dxi + N j dt) . 
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It turns out to be technically simplest to work in the gauge: 

6<p = cp, hij = e 2<T ((5y + 7y) , dijij = , j lt = . (10) 

This gauge is different from the unitary gauge 5<p — used in [7]. However, we will comment on their relation in due 
course. 

To relate this gauge choice to the physical (conserved) curvature perturbation, we can use the 5N formalism [l7j |. 
In this formalism, one can relate the gauge invariant curvature perturbation £ to the number of e-foldings since the 
time of horizon crossing 

C(t, x) = N(t, x) - N (t) = SN , 

where 

N(t,x) = J H(t',x)dt' , (11) 

and No denotes the background value without the perturbation. We can then view N as a function of the scalar field 
perturbation evaluated at the time i*. We can then write, 

:Kr dN ld 2 N 2 

m = W^ + 2W/* + - (12) 

The first derivative follows directly from (fTTjl : 

ON _ _#* 

d<l>* <j>* 

where, as usual, the star denotes evaluation at t*. 

Higher derivatives of N will be suppressed by powers of the slow roll parameters. For example, by using the 
equations of motion ^ and ([5]), and in the limit where e = 0, one can show that 



It is then easy to show that the extra terms in the expansion (|12[) will give corrections to //vl in powers of ~ 
(slow roll) /Cg . We assume that the slow roll parameters are small such that (slow roll) /c 2 <C 1 ■ In any case, we will 
be working in a de-Sitter background where these corrections are exactly zero. 
The two point function of the curvature perturbation is given by, 

<Cjf(*)C*(t)> « ^(n(t*)^(t*)) = (2nfS^(k + k')P < (k) . (13) 
Similarly, the three point function can be written as 

= (2vr) 3 ,5^(fci + k 2 + k 2 ) f-px^Y^PdWikj) . (14) 

\ ' i<j 

In the last line we have defined the non-gaussianity parameter Jnl, following the conventions in 18]. 
To linearized order, the solution of the constraint equations can be parametrized by two scalars, 

N = 1 + SN , N { — diX ■ 

After many integrations by parts, the quadratic Lagrangian for the Einstein-Hilbert, scalar and Gauss-Bonnet part 
of the action read respectively, 

4h = \ J e 3a H5N (3HSN + 2V 2 X ) , (15) 

Sf ] = 1 J e 3a [dlPip 2 + 2X 2 8 X P + 2SN {d^Ptp + Xxd x P) + X x {2d <t> d x Pv + X x d x P)] , (16) 
S^l = - J e a 8V(H 2 8N\7 2 ip + J e 3a {-48V!H 3 <j)5N 2 + 12(H 2 + H)H 2 V{ 

f(V[H - 4>V{') - V{ip +2AH 2 5N tpV"H<j) + V{{-<t>V 2 x + Hip) X . (17) 



-8i? 2 V 2 x 
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where 



X 1 = (f>{ip - <j>5N) , (18) 
'' [ 3<f) 2 8N 2 - e - 2<T (V^) 2 - 2<pdiipdiX ~ tySNip + ip 2 } . (19) 



X 
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No approximations have been made in deriving Eqs. (|15p - (|19[) . 

As pointed out in the previous section, in the slow roll limit we can ignore the terms involving the second derivative 
of the potential V\. Varying Eqs. (fT"5|) - (TTTj) with respect to the constraints we obtain, 



c s g <p 
2a(l + §<?) H 



5N « _ / , a 3 ^ ~~ , (20) 



3 ^ 1 



ffCs :e- 2CT VV+^(^-#iV) 



(21) 



2aff r 2ac s 

where we have kept only the leading terms in the limit c s — > [22] . Therefore, we see that in this limit, 

(57V ~c s , x~~- (22) 

This order-of-magnitude estimate is important to determine which terms in the action survive the small speed of 
sound limit. 

Inserting the solution for the Lagrange multipliers (|2"D|) and (|2"Tj) back into the quadratic actions (|15[) - (fT7|l and 
taking the small speed of sound limit, we obtain the quadratic action for the scalar fluctuations: 

krn = J e 3 V(.9) ~ ^(V^) 2 e - 2 1 (23) 

where 

25(1 + gf 



f(g) 



a 2 (2 + 3g) 2 ' 



. 2 (l + 2. 9 )(2 + 3ff) 2 
Cs 2(1 + .g) 2 ° s ■ 

Note that the quadratic action for the fluctuations is non-trivial even in the de-Sitter background that we are consid- 
ering. Moreover, we see that g must be positive in this limit in order to give the correct sign for the kinetic term in 
Eq. HMD- 

The dispersion relation for the quadratic fluctuations is precisely of the relativistic type u> ~ k. Such fluctuation 
spectrum around de-Sitter was first found in |7|] by studying the effective action for the scalar fluctuations directly. 
We then see that our approach gives a physical interpretation to such fluctuations: they are generated by the Gauss- 
Bonnet coupling of the inflaton. The precise dictionary between our variables and those of Ref. [3] is the following: 
( = —Hit = (H/(j))tp, where it is the "Goldstone Boson" of [7(|. Moreover, in our gauge, the extrinsic curvature and 
the Lapse can be written as 

SKi^-didjX, 5N^-±(l + g m ) 

It is then easy to see that the terms that give the relativistic dispersion for the quadratic fluctuations come precisely 
from couplings of the form (1 + g 00 )5Kl in the effective action. These were the terms studied in [7j]. 
The quantization of the perturbations using the action f|23[) proceeds in the standard way. We write, 



f(j,x)= / 7^^(r) 



e 



ik-x 



(2tt) 3 

where the Fourier transform is written in terms of the standard harmonic oscillator operators 

^( r ) = u k( T ) a l + u H T ) a -k ■ 



Here we are using conformal time defined by dt = e CT c?r. For de-Sitter space, e CT = 



-l/(Hr) where r S (-oo,0]. 
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The properly normalized Bunch-Davies vacuum is given by [191 ], 

u r = (1 - ikc s r)e lk ^ T . 



The power spectrum follows from the definition (|13p : 

4- 1 1 

(24) 



_ ff 2 (2 + 3g) 2 1 _^ 1_ 
CW " 8 ff (l + «?) 2 c s fc 3 " C fc 3 



Note that, in the slow roll limit, we have a scale invariant spectrum. Therefore, observations of the scalar tilt do not 
put constraints on the value of the Gauss-Bonnet coupling. 

III. NON-GAUSSIANITIES 

To study the non-gaussianities for this model, we need to expand the action to cubic order. In this section we will 
show only the leading terms in the limit of a small speed of sound. Moreover, we will always ignore the mixing with 
gravity. To estimate the size of the various terms, the limits (|2"2"|) are useful. It turns out that the cubic action is of 
order ~ l/c s . There are many integrations by parts necessary to put the results in a simple form. The final result for 
the cubic couplings takes the form, 

HmS® = -\J ^SNldidjxdidjx-^xf] , (25) 

(3) f ( „o „ 1 



Hmfi^ = J e^\X 1 X 2 dj c P+-Xtd 6 x P\ , (26) 
c lrm4 3 B - - / e '^ SN [M^9 3X -(VV)(V 2 x)] - | f e 3a gSN [d^xd^x ~ (V 2 *) 2 ] 



9 -^ dlV d lX V\ + \ J e 3 ^ [WxMiX ~ (V 2 X) 2 ] , (27) 



where, 



Xi = [ip 5N 



lim X 2 = -I e - 2CT (V^) 2 - —dupdiX ■ 

c s ^0 2 C s 

In the usual k- inflationary scenario (g — 0), the Einstein-Hilbert action Eq. (f2"5|) gives a contribution to the tri- 
spectrum which is of higher order in slow roll: /nl ~ e2 / c s- Therefore, this term has been ignored in previous 
calculations. However, with a non-zero Gauss-Bonnet coupling, we have seen that fluctuations can exist even in a 
de-Sitter background. Therefore one can have j»e. In this case, which is the main focus of this paper, all terms in 
Eqs. (|23|) - (|2T|) are equally important. 

In the language of the effective field theory of Q, we can readily identify some of the cubic couplings. For example, 
the Einstein-Hilbert term can be written as S^h ~ SO- + - (<^Q) 2 )- However, since our gauge choice is 

different from [7j, not all terms in (|25[) - (|2T[) can be written only in terms of the Lapse and the extrinsic curvature. 

In order to calculate the scalar three point function we follow the standard procedure. To leading order in the cubic 
perturbation, the three point function is given by 



J —OC 

« -if dr , e^')(0|[^ i (0)^(0)^(0),ff in t(r')]|0) > 



where in the last line we have changed the integration to conformal time, and made the usual late time approximation. 
The interaction Hamiltonian is given in terms of the cubic Lagranian density: H mt = —£^ 3 '- 
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The cubic Lagrangian density can be written in Fouriers space as 

d 3 k x d 3 k 2 d 3 k 3 



(2n) 3 6 ( - 3 Hk 1 + k 2 + £ 3 )£ (3) (fci, k 2 , k 3 -t) 



(2tt) 3 (2tt) 3 (2tt) 3 

Then, it is not hard to see that the three point function can be written as, 

( ni (t) n2 (t) ni (t)) = ^(2n) 3 5^(k 1 + k 2 + k 3 )f[\u^(0)\ 2 f ^-& 3 Hh,k 2) k 3 ;x) 

J -co aX 



c.c. + perms. 



where we have defined the integration variable x = c s t. Moreover, C^(ki, k 2 , k 3 ; x) is calculated by taking the cubic 
Lagrangian density that follows from Eqs. P5|) - (f2"T)) and replacing the fields <p using the rules: 



2 2ikx 



<p->(l- ikx)e lfcx , ip -> -Hk 2 x z e 

where k = |fe|. Moreover, all spatial derivatives are replaced by their fourier transform: di — > i{k)i. Using the 
definition of /tvl in Eq. (] 14[) . we can write 



fNL(h,k 2 ,h) 



5a 3 V c 

3 

i '"i 



6c 3 E^ 



. I ^,£(3)/^ ^ 2 ^ 3;a ,\ _|_ c c + perms. 

J-oo ^ 



(28) 



In doing the integrals, one projects to the correct interacting ground state by rotating slightly the contour so that 
t — > (1 — ie)r. With these considerations, it is straightforward to calculate the three point function. As an example, 
let us work out the Einstein Hilbert term, Eq. (|25[) . The Fourier transform Lagrangian density can be written as, 



where, 



and 



^EH( fc i> fc 2,/c 3 ;a;) 



SN(k) 



Hx J 



(k 2 ■ h) 2 - k 2 k 2 57V(fci)x(fc 2 )x(fc 3 )e^ 



■^k 2 x 2 , x(k) = ^C\-2 + b-2a 2 (l-ikx)} 
2a Aaac* 



c s =ac s , K = k x +k 2 + k 3 
Therefore, the time integration gives, 



a 2 = 



2(1 + gf 



(l + 2.g)(2 + 3.g) 



6 = 



1 + 1.9 



J — oo 



ib 3 c 2 s k 2 
32aH 2 K 5 a 3 



(fca • h) 2 - k 2 2 k\ 



A8k 2 k 3 a 4 + 6 (2a 2 -6 + 2) K(k 2 + k 3 )a 2 + (-2a 2 + 6-2) K 2 
Using the definition (128|) we find the contribution to Jn l of the Einstein-Hilbert action: 

56 2 (36 - 2)kl \{k 2 ■ k 3 ) 2 - k\k 2 3 



f (EH) 
J NL 



{6 [(2a 2 -6 + 2) Kk 3 



96a 4 (6-2) 2 £ 2 K 5 kf 
+k 2 (8k 3 a 2 + (2a 2 -6 + 2) K)] a 2 + (-2a 2 + 6 - 2) 2 K 2 } + perms. 



(29) 



In calculating the contribution from the scalar sector, one encounters the question of how big is d x P compared to 
d x P (note that both X\X 2 and X 3 scale as ~ 1/c 3 ). There is no way of knowing this without a detailed form of P. 
We will simply introduce a new variable A and write, 



x ^ a ^d 3 x P_d x Pd 3 x P 



d%p (d 2 x py 



(30) 
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Then, following the same steps as in the Einstein-Hilbert term, it is easy to show that the contribution from the scalar 
sector to f^L is given by 



M) _ 5(2 - b)\k 2 k 2 k 2 



f w _ 
J NL ~ 



5(fc 2 ■ k 3 )k 2 
+ 12a 4 (2- 6)c 2 if 3 
+perms. 



[2(6 - 1) (K + 2k 2 ) k 3 a 2 + (2a 2 (b - 1) - (6 - 2)6) K (K + k 2 )] 



(31) 



Finally, the contribution from the Gauss-Bonnet term gives, 
56 2 



,(GB) 
J NL 



| (b - 2) ((k 2 ■ k 3 ) 2 - kjkfj [I8k 2 (8k 3 a 2 + (2a 2 -b + 2)K)a 2 
2) K)} k 2 



48a 4 (6-2) 2 c ;2 if 5 j-.kf 
+ (6a 2 - 36 + 2) K (6k 3 a 2 + (2a 2 - b 

+2K(k 1 ■ k 2 )k 2 k Y U ((2a 2 -6 + 2) Kk 3 + k 2 (6k 3 a 2 + (2a 2 -6 + 2) K)) 



+2K 2 {k 1 ■ k 2 )k 2 2 ((2a 2 
+perms. 

The total contribution to /jvl is then 



6 + 2) Kk 3 + k 2 (Ak 3 a 2 + (2c 



2)K)) 



-2a 2 
-2a 2 - 



6-2) K" 



h -2) 2 if 2 ]} 



, _ , (EH) M ,(GB) 
JNL - J NL + J NL + J a 



I NL 



(32) 



(33) 



where the different terms are given in Eqs. (|29|) . (f3Tj) and (|32|) . Note that the wave- vector dependence of Jml can be 
written in terms of the magnitudes ki by using, 



ki ' k*j 



1 



(k 



kj) , 



J ~ 2 ■" 

If we take the limit of vanishing Gauss-Bonnet coupling g — > 0, only the scalar contribution survive and we get, 



lim f NL 



5 A/b-^ k^ 
W^E~kf 

5 A/c^ k<2 k^ 



24^ E . fc 3 
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JsT (K + fc 3 ) + ^2 (K + 2fc 3 )] + perms. 



2c 2 if3 



2if 2 



(34) 



One can check that Eq. (|34|) coincides with the Jnl calculated in Q in the case of k-inflation, and in the slow 
roll and small speed of sound limits (28[. If one is interested in a finite speed of sound, one then simply replaces 
1/c 2 — > (1/c 2 — 1) in the second term of Eq. (f3"4"| . 

So far we have assumed that the speed of sound is small. However, another way of getting large non-gaussianities 
is to make A large. In this case, /nl will be dominated by the first term in Eq. (|31[) : 



5(2-6)Afc 2 /c 2 fc 2 



(35) 



Note that the dependence on the speed of sound is exact in this case (see Q and footnote [28]). Therefore, in this 
general class of models, we can parameterise the leading contribution to f^L by three parameters: c s , A and g. In 
the next sections we will see how observations put constraints on these parameters. 



A. The Shape of Non-Gaussianities 



Studying the shape of fjq l amounts to calculate how it depends on the size and shape of the a triangle with sides 
fcj. To get an easier visualization, we will follow [l8l |. and plot Jnl as a function of the two ratios: x 2 = k 2 /k\ and 
x 3 — k 3 /ki. Note that the dependence of k\ drops out. We will begin by studying the effects of the Gauss-Bonnet 
coupling. Therefore we set the extra parameter A = for the moment. 
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1.0 xi 




x 3 



FIG. 1: The non-gaussianity parameter /jvl as a function of X2 = fe/fti and X3 = fe/fei for p = A = 0. We have normalized 
Jnl so that it is equal to one in the equilateral limit ki = kz — k$. 

The shape of Jnl for A = g = is shown in figure 1. We normalize /jvl so that it is equal to one in the equilateral 
limit k\ — hi — k^. However, we respect the sign that follows from the definition given in the previous section. 
Moreover, in order to avoid overcounting configurations, we set fjsrL = if xi and X3 do not obey the inequalities, 
1 — X2 < £3 < X2 ■ The shape shown in figure 1 is what we expect in usual k- inflation (with A = 0) . 

If we turn on the Gauss-Bonnet coupling, we get a slight deformation of this shape. However, it is easy to show that 
/jvl is bounded even if we take g — > 00. The deformation of /jvl is clearest if we plot the change of the normalized 
non-gaussianity, 

Af NL = f%>r( g )-f%>r(g = o) , 

where fjlf™ means that we have normalized /ml as described in the previous paragraph. In figure 2 we plot A/jvl 
for g = 3. In the next section we will see that this value of g is within the current observational limits. It is 




FIG. 2: Deformation of Jnl due to the Gauss-Bonnet coupling for g — 3 (A = 0). 

interesting that, for all values of g the deformation of Jnl obeys the bound |A/jvl| 10 _2 |/7VL(ff = 0)| ~ 10 -2 /c 2 s . 
However, this is about the same order of magnitude as the contributions coming from the slow roll parameters which 
are generically of order ~ e/c 2 . ~ \n s — l|/c 2 ~ 10 _2 /c 2 , where n s is the scalar tilt. In the most optimistic scenario, 
where 1/c 2 ~ 10 2 , these deformations will represent a change of order |A/jvx| ~ This is about the lower end 

of the detectability threshold of non-Gaussianities 20]. Therefore, we conclude that the Gauss-Bonnet coupling does 
not produce a measurable deformation of the shape of Jml from that of k-inflation. 

Let us now turn our attention to the effect of the parameter A. To estimate the biggest possible contribution from 
this term, we assume A ^> 1. Then, we plot the difference 

Af NL = f n N T\g = 0, A - 00) - f%>r(g = 0, A = 0) , 
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FIG. 3: Deformation of /ml due to the Gauss-Bonnet coupling for g — > oo (A = 0). 



where 



jTi m (<? = 0,A^oo) = 



The resulting plot is shown in figure 4. We can see that the deformation due to A is an order of magnitude larger 




1.0 



0.8 0-6 



0.4 0-2 



FIG. 4: Deformation of Jnl for A > 1. We have set g = in this plot. 



than the rest of the deformations. Therefore, it is potentially observable. Although the plot of figure 3 is for A — * oo, 
we have found that the deformation is of order |A/$£ ro | - ©(KT 1 ) for |A| > |1 - c^lO 3 . 

If a large Jnl is observed in the CMB, one could test whether the dominant effect is a small speed of sound or 
large A, by looking at the deformation around the basic shape with g = A = 0. However, such test will not put any 
constraints on g since we have seen that its effect on the shape of Jnl is negligible. Nevertheless, there is another 
effect of the Gauss-Bonnet coupling that we have not studied so far. In the next section we show that the Gauss- 
Bonnet coupling amplifies the spectrum of gravitational waves. We will see that this enhancement allows to take a 
smaller speed of sound, and so in an indirect way, the Gauss-Bonnet coupling can also lead to an amplification of the 
non-gaussianities. We will show that current WMAP data can already put constraints on the value of g. 



IV. THE GRAVITATIONAL WAVE SPECTRUM 



In this section we will study the gravitational wave spectrum of the Gauss-Bonnet coupling in the slow roll limit. 
We will only consider the quadratic fluctuations. Similar studies have been done in [2l| for a field with canonical 
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kinetic term. Note that the axion coupling in the action ([T]) does affect the gravitational wave spectrum. However, 
the effects of this term were understood in Q and (l4| . The axion coupling introduces a helicity dependence in the 
tensor power spectrum. We will ignore this term in the following. 

The Gauss-Bonnet term, however, has a very different effect. It allows for both a large amplitude of gravitational 
waves and large non-gaussianity. This can be traced back to the fact that g can be much larger than the slow roll 
parameters. Moreover, as we will show below, a large g with a fixed tensor/scalar ratio, also leads to a smaller value 
of the speed of sound. 

The quadratic actions for gravitational waves are given by, 



Q7 
°EH 



3 a 



2{3H 2 + 2H)% 



-2a 



SJ = - 



e 3CT P(X,0) 7 2 , 



C7 
°GB 



l J e ' CT [9il - ^H 2 ( Vl + g(2 -2e + OH 



-2a 



0?i+6?)(V7ij) 2 



In the slow roll limit, the total quadratic action simplifies to give 

1 



S 7 = 



e 3CT [{l+g)%-e- 2 °{V ll3 ) 2 



where we have used the equations of motion ((4]) and |5]). 
The power spectrum is then easily calculated: 



where s, s' labels the two helicities of the graviton. Using the scalar power spectrum Eq. (f24|) . the tensor/scalar ratio 
takes the form, 



' 7 



lQc s g{l + gf' 2 
(2 + 3 5 ) 2 



(36) 



In the limit where the Gauss-Bonnet coupling vanishes, the dominant contribution to r will come from the slow roll 
parameter e [j| 



lim r = 16c s e . 



(37) 



Therefore, we see that a small speed of sound suppresses the gravitational wave amplitude in this case. However, in 
the case of a non-vanishing Gauss-Bonnet coupling, we see from (l36|) that g can be relatively large to compensate 
for a small speed of sound. In this way one can produce a large amplitude of gravitational waves and a large value 
of /jvl- Of course, as we discussed in the previous section, one can also produce large non-gaussianities with a large 
value of A. In this case, the speed of sound does not need to be large. 

We now want to know if we can constraint the parameters c s , g and A using current CMB data. The current 
limit on the tensor/scalar ratio is r < 0.20 Q. Therefore, from Eq. p6|) we see that g cannot be too large without 
decreasing the value of the speed of sound. The allowed parameter space for c s and g is shown in figure 5. 

Note that for non-trivial values of the Gauss-Bonnet coupling (g > 0.1) we need a small speed of sound. This will, 
in turn translate to large non-gaussianities. Moreover, a small value of g will be very hard to disentangle from the 
contribution of the slow roll parameters. Therefore, we see that the most interesting region is that of a small speed 
of sound, and hence large Jnl- This is consistent with the approximations made in the calculation of Jnl- 

We now want to constraint A and g combining limits on r and on Jnl- To do this, we eliminate c s in terms of r 
and g using Eq. (|36[) . We can then insert this value into /atl, Eq. (|33|) . and write the equilateral /jvl as a function 
of r, g and A. In the equilateral limit, the formula for /jvl simplifies to, 



foq 
J NL 



?2 



5A(1 + 2 5 ) 2 (2 + 3 5 ) 5 (2ei + 1) (ei (ei (ei (72ei + 379) + 690) + 516) + 136) 



486(1 + g) 3 



1296 (ei + l) 4 (3ei + 2) 



(38) 



where, 



r(2 + 3g) 2 
16 5 (1 + g) 5 / 2 
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FIG. 5: Observational constraints on g and c a . The shaded region has been ruled out by the observations. 
3.0 
2.5 
2.0 
$ 1.5 
1.0 
05 
0.0 




•100 



-50 





A 



50 



100 



FIG. 6: Observational constraints on g and A. The shaded region has been ruled out by the observations. 



< 



250 Given these constraints, we 



Current WMAP bounds on the equilateral non-gaussianities are roughly |/^ 
show the allowed parameter space for g and A in figure 6. 

One can see a degeneracy point around A ~ 2 where one can take large values of g. This is due to the fact that 
we can have cancelations between the two terms in brackets in Eq. (|38|) . However, this is a very fine tuned situation. 
Nevetheless, we see that most of the parameter space for g is quite well constrained by current observations. This is 
specially true for models with large |A|. 

The bound on A is, however, not very good. The absolute limit comes from taking g — > and using the limit of r 
given in Eq. (|57|) . The equilateral Jnl simplifies to 



lim „ Inl = 



_ 320e 2 (4A-51) 



243r 2 



(39) 



Since e is enters linearly in the scalar tilt, one has that e < 10 2 given current WMAP observations. Using the limit 
on equilateral /jvl and r given above, one gets a rough bound 



|A| < 10 4 



(40) 



So far we have concentrated in giving bounds on g and A. However, one can ask if it is possible to break the 
degeneracy between these two parameters and put more precise constraints on a non-zero Gauss-Bonnet coupling. 
The most optimistic scenario would be to detect a large /jv l an d the deformation due to A described in the previous 
section. In this case one would confirm that |A| > |1 — c 2 |10 3 . This, combined with limits on r (see figs. 5 and 6) will 
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put tight constraints on g. Nevertheless, in most models A is not an independent parameter. In this case we can put 
more precise constraints on the Gauss-Bonnet coupling. We illustrate this point in the next section for the case of 
DBI inflation. 



V. AN EXAMPLE: DBI INFLATION 



So far we have discussing the constraints on the Gauss-Bonnet coupling in a model independent way. However, in 
order to break the degeneracy between the parameters g and A we need to consider particular models. In this section 
we will study DBI inflation [22j | . This is a string-inspired model where the inflaton encodes the position of a D-brane 
on a warped compactification. In this model, the function P{X, <fi) takes the form 



p(x, 0) = -ftf)- 1 y/i - 2XJW) + f{4>T l - v{4>) , 

where /(</>) is related to the warp factor of the compactification. 
Then, it is easy to show from the definition of A, Eq. (|30|). that 

Adbi = 3 . 

In this model we do not have any degeneracy between the parameters and we can constraint g directly. In figure 7 
we plot the the absolute value of the equilateral Jnl as a function of g and the tensor/scalar ratio r. We see that in 
DBI inflation, the Gauss-Bonnet coupling is constrained as 



Sdbi < 3 . (41) 



Moreover, we see that there is a very interesting region of parameter space where we can have both large /jvl and a 
large amplitude of gravitational waves. 




0.20 



FIG. 7: Equilateral Jnl for DBI inflation, as a function of the tensor/scalar ratio r and the Gauss-Bonnet coupling strength 
g. The observationally allowed region is for r < 0.20 and \f^ L \ < 250. 

So far we have assumed that the slow roll parameters are very small compared to any other scale in the problem. 
This will translate in some constraints on the potentials V, V\ and /. In the slow roll limit, the equations of motion 
0} and ||SJ) reduce to, 

§2 = W , (42) 
--l]+V(4>) = 3ff 2 (l+.g), (43) 
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where, 

c s = y/l - 2Xf . (44) 



To solve these equations, one first solves for 4> from Eq. (|42|) and substitute it on Eq. (|43|) . This gives a Hubble 
parameter that is a function of the scalar field H = H ((/)). We then substitute this back in (14"2")) to solve for <j> as a 
function of time. We will not do this in details, since we are only interested in deriving general conditions on the 
potentials. 

In the limit of a small speed of sound, one has cp 2 ss 1/(2/). One can then show that 

c s U) pb == (47) 

From Eqs. (f4"5)) and (jirJ)) . we see that -H" obeys a cubic equation. However, we do not need to solve this equation 
explicitly to derive the slow roll conditions. Using Eqs. (|45l) - (|47| along with {8} and ([9]), one can derive the following 
relations between the slow roll parameters: 

. » H^) * , 9 9 



2V{(j))H 3 + 2gHg 

vi' it) i> 



Hg ' Vi{4>) H 



(49) 



|S + l) + 3, (51) 



Hc s H\V{(<P) 2/(0)7 

Our results have assumed that all slow roll parameters are much smaller than 4>/H. Therefore, to be consistent we 
should impose the conditions: 

^«1 ™«1 ^M«l ^M<1 (52) 

m « ' « ' ^ ' u(^) « ' (52) 

where the last condition follows from e/(He) <C 1. 

For an AdS throat, the warp factor is /(</>) = A/(/> 4 , where A is the 't-Hooft coupling of the dual gauge theory, and 
A ^> 1 [22j]. Therefore, the first condition in (|52")l gives, ^ > 1, In other words, the D3-brane must be in the UV 
region of the AdS throat. Note that this is still consistent with the bounds studied in [22j, so that one can ignore the 
backreaction on the geometry. 

If we want to be in the interesting window of large /nl and observable gravitational waves, we need g ~ 0(1). In 
general, this condition will require a very steep potential V\. In fact, using Eq. (|36[) for the tensor/scalar ratio, along 
with Eq. (|4"2"|) for <j>, one can show that in this regime V{{<ft) ~ 10 10 . Since, for the AdS throat, one is interested in 
^ > 1, one could perhaps realize such large value of V[ ((/>) with a power-law potential V\{4>) ~ (j) n . Note that the 
slow roll conditions in (|52")l will automatically be satisfied for large values of cf>. Wether such scenario can be realized 
in a controlled string theory construction, is beyond the scope of this paper. 



A. Comments on the Lyth Bound 

We have seen that a Gauss-Bonnet coupling enhances the amplitude of gravitational waves. It was shown in [lj| 
that, in the context of slow roll inflation, an observable amplitude of gravitational waves would require a ultra- 
planckian displacement of the inflaton. This is known as the Lyth bound. Quite generally, we can write the bound 
as 



Ac/) > \<j>/H\AAT, 



(53) 
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where AAT is the number of e-foldings in which the scales of interest in the CMB today exit their horizon. This is 
usually taken as AAf ~ 4.6 15]. Using Eqs. (j4"2)) and (f3"6"| in ([55]) . one can write the Lyth bound in terms of the 
tensor/scalar ratio and the Gauss-Bonnet coupling: 



Acj> > 0.36183 



'Sg + 2 



(.9 + l) 3/4 (ff(6 5 + 7) + 2)i/4V 0.07 



0.693566 



.9 



1/4 



0.07 



(54) 



0.07 is considered to be the lowest limit of 



where in the last step we have shown the limit of large g. Moreover, 
detectability for gravitational waves . 

We can see from Eq. ([54"]) that one can, in principle, violate the Lyth bound if we have a large Gauss-Bonnet 
coupling. However, a large value of g, with fixed r requires a very small speed of sound (see Eq. (j36j) ). We have seen 
that current WMAP limits on Jn l already constraint the value of the Gauss-Bonnet coupling to Eq. (|4T|) . With this 
value we get A(f> > 0.5, which is roughly the same bound as in slow roll inflation. Nevertheless, if we are working 
with an AdS throat, we need <j> ^> 1 and so A0/0 < 1. So we see that the fractional change in the scalar field is very 
small. This means that, in order to make CMB predictions within this class of models, we only need to know the 
potentials in the large <p limit. It would be interesting to see wether one can violate the bound in some other model 
with a Gauss-Bonnet coupling which allows ~ 1. 



VI. CONCLUSION 



In this article we have studied WMAP constraints on modifications of gravity due to a Gauss-Bonnet coupling in 
single-field inflation. This is the most general modification of gravity that leads to second order equations of motion, 
and that also affects the spectrum of scalar fluctuations. We showed that in the slow roll limit, and for a very general 
class of models with action ([T]), a large Jnl can be written in terms of three parameters: c s , A and g the Gauss-Bonnet 
coupling. 

We found that the Gauss-Bonnet term has little effect on the shape of non-gaussianities, but it amplifies the tensor 
power spectrum. Thus, given current limits on the tensor/scalar ratio r, we found that large values of the Gauss- 
Bonnet coupling would require small values of the speed of sound, and hence large non-gaussianities. Using current 
WMAP limits on r and /nl we were able to constraint the parameter space of such models. 

To give better constraints on g we studied a particular model: DBI inflation. In this case we obtained a precise 
bound on this coupling, Eq. (|4"Tj) . Moreover, we saw that in this model, a non-zero Gauss-Bonnet couplings leads to 
an interesting observational window with both large non-gaussianities and a large amplitude of gravitational waves. 
We also studied how the conditions for the smallness of the slow roll parameters translate to constraints on the scalar 
potentials. Possible violations of the Lyth bound were also studied. We found that the bound can be violated for 
large values of g. However, for DBI inflation one has an observational restriction of g < 3, and so the bound is 
roughly A(j) > 0.5. Nevertheless, we found that for an AdS throat the fractional change in the scalar field is very 
small: Atfi/cj) <C 1. It would be interesting to study the higher curvature corrections to the DBI action, to see if one 
can realize a Gauss-Bonnet driven inflation in a controlled way. Some of these corrections were derived in [23| . 

Another interesting aspect of the Gauss-Bonnet term is that scalar fluctuations can exist in a de-Sitter background. 
Moreover, they have a relativistic dispersion relation, unlike Ghost Inflation [161 ]. This matches with the new de-Sitter 
limit found in 7] using the effective action of the scalar fluctuations. One might wonder if there are other de-Sitter 
limits of inflation. In 7}, it was argued that the answer is negative as any such limits will not make sense as an 
effective field theory for the fluctuations. However, the authors of [7[ only considered models with one scalar degree 
of freedom. If we ask about other modifications of gravity, it is well known that one needs to add more degrees 
of freedom to the theory. It would be interesting to put constraints in other types of modified gravity using large 
non-gaussianities . 
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